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In this paper we compute the entanglement, as quantified by negativity, between two blocks of 
length La and Lb, separated by L sites in the one dimensional spin-1 AKLT model. We took the 
model with two diflerent boundary conditions. We consider the case of A'^ spins 1 in the bulk and one 
spin 1/2 at each boundary which constitute an unique ground state, and the case of just spins 1, even 
at the end of the chain, where the degeneracy of the ground state is four. In both scenarios we made a 
partition consisting of two blocks A and B, containing La and Lb sites respectively. The separation 
of these two blocks is L. In both cases we explicitly obtain the reduced density matrix of the blocks 
A and B. We prove that the negativity in the first case vanishes identically for L > 1 while in the 
second scenario it may approach a constant value N = 1/2 for each degenerate eigenstate depending 
on the way one constructs these eigenstates. However, as there is some freedom in constructing these 
eigenstates, vanishing entanglement is also possible in the latter case. Additionally, we also compute 
the entanglement between non-complementary blocks in the case of periodic boundary conditions 
for the spin-1 AKLT model for which there is a unique ground state. Even in this case, we find that 
the negativity of separated blocks of spins is zero. 



During the last years, in the disciplines of quantum in- 
formation theory, quantum many-body physics and sta- 
tistical mechanics, the study of purely quantum effects 
like entanglement has become important. The idea of 
entanglement was first introduced by Schrodinger Ij in 
the early days of quantum mechanics. Entanglement is 
a phenomena where two (or more) quantum systems are 
linked together and their description cannot be done sep- 
arately, disregarding their spatial separation. 

Entanglement plays a fundamental role in some quan- 
tum mechanical systems, being important in the descrip- 
tion of quantum phase transitions, [H, topological order 
[31 13] and even macroscopic properties of solids [5] . For 
pure systems, the measure of bipartite entanglement is 
given by the entanglement entropy, or von Neumann en- 
tropy, defined as S[p\ — Trplnp, where the state is char- 
acterized through it's density matrix p. 

Unfortunately for mixed states, the von Neumann en- 
tropy is not an appropriate measure of quantum entan- 
glement. One is thereby forced to use other measures 
of entanglement for mixed states [3 [7j. In this letter 
we focus on a particular measure of entanglement called 
the negativity, which stems from the Peres separability 
condition 0, and first proven to be a bona fide quan- 
tifier of entanglement in [3]. If one is to quantify the 
genuine quantum correlations or entanglement between 
non- complementary parts of a many body system (as op- 
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posed to mutual information, which only quantifies total 
correlations) , one has to resort to a measure such as nega- 
tivity. Unfortunately, such calculations have been found 
to be only numerically tractable so far [10], except for 
infinite range models [TT] . 

Given a system E, we define two subsystems A and B, 
characterized by the density matrix Trcp = pas, where 
we have traced away the degrees of freedom C, who lie 
outside the A and B subsystems. Negativity is defined 
then as the sum of negatives eigenvalues of the partial 
transposed density matrix p^g, where the transposition 
is done in the A subsystem. The existence of nonzero 
negativity is signal of entanglement, but, as pointed out 
in [13^ , vanishing negativity does not imply zero entan- 
glement. 

In this paper, we present the result for the negativity 
of two blocks defined as subsystems of the AKLT ground 
state. In quantum information, entanglement between 
separated systems is a valuable resource. In this sense it 
is important to calculate the entanglement between sep- 
arated (i.e. non-complementary) blocks of a system. It is 
therefore important to examine whether such an entan- 
glement resource is present in the AKLT model whose 
ground state has already been demonstrated to be an 
ideal resource for measurement based quantum wires [TH] 
in ID, and quantum computation |19] in 2D. It is a re- 
alistic model of spin systems in the sense of being short 
range and models the gapped nature of integer spin mod- 
els. Perhaps its most attractive feature is its analytic 
solvability and thereby it holds the promise of the first 
ever analytic computation of negativity for a short range 
model. 
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AKLT model 



We will study two different scenarios for the problem. 
In the first case, we will compute the negativity in the 
AKLT model for a bulk of spins 1 and one spin 1/2 at 
each boundary. This system has a unique ground state, 
as was shown in jl6j for the general case. In the second 
scenario, we study the case of a bulk and the boundary 
made up of spins 1. This system should be easier to real- 
ize experimentally than the previous one, but it has some 
subtleties as it's ground state is four fold degenerate. 



S 



albi 



a-jb] 



±1 _ _i_c±x rc+ 



with \si,sf\ = ±stb,,, \st.s-\ 



-bib,), 



(5) 



where a and b are two sets of bosonic creation operators. 



with the usual commutation relations [a;, Oj] — [bi,b 
6ij, [ai,aj] — [bi,bj] — and correspondingly for and 
b^ . This two sets commute in each and every lattice site, 
i.e. [ai,bj] = [al,bj] = 0. In terms of these variables, the 
ground state can be written as 



A. Spin 1/2 at the boundary 

The one dimensional AKLT model [15] that we will 
consider consists of a chain of N spin-Is in the bulk, and 
two spin-1/2 on the boundary. The location where the 
spins sit are called sites. We shall denote by Sk the vector 
of spin-1 operators and by Sb spin— 1/2 operators, where 
b = 0,N + 1. The Hamiltonian is: 



N 



H — i?Bulk + Ho,! + I\-N,N+1, 



(1) 



where the Hamiltonian corresponding to the bulk is given 

by 



N-l 



-f^Bulk = PjSj + Sj+i). 



(2) 



i=l 



Here P{Si + Sij^i) is a projector onto spin 2 states, given 
by 

P{Su + Sk+i) = \ {2.Sk ■ Sk+i + {Sk ■ Sk+if + 2) , (3) 

and the sum runs over the lattice sites. The boundary 
terms H describe interaction of a spin 1/2 and spin 1. 
Each term is a projector on a state with spin 3/2: 



2 ^ 2 ^ 

IIo,! = -(l + So-5'i), Hjv,Ar+l = -(I + S'at-sV+i). (4) 

In order to construct the ground state |VBS) of ([T]) we 
can associate two spin 1/2 variables at each lattice site 
and create the spin 1 state symmetrizing them. To pre- 
vent the formation of spin 2, we antisymmetrize states 
between different neighbor lattice sites. Doing this we 
are sure that this configuration is actually an eigenstate 
of the Hamiltonian ([T]), with eigenvalue (i.e. the pro- 
jection of |VBS) on the subspace of spin 2-states is zero). 
Noting that the Hamiltonian ^ is positive definite, then 
we know that this is the ground state. It is possible to 
write down a compact expression for this VBS state us- 
ing bosonic variables. Following [12], we make use of 
the Schwinger boson representation for SU{2) algebra at 
each site j, namely: 



|VBS)^n(°M+i-«li^l)|0)- 



(6) 



i=0 



where |0) = 0^^^^^ \Oa,j) «) |Ofc, j). The state \Oa,j) is 
defined by aj\Oa,j) = 0, and it's called the vacuum state 
for the set of operators a. \Ob,j) is defined similarly for 
the set b. In [16] the authors prove that this ground state 
is unique for the Hamiltonian (IT]), then we can construct 
the density matrix of the (pure) ground state 



P = 



|VBS)(VBS| 
(VBS|VBS) ■ 



(7) 



This is a one dimensional projector on the ground state 
of the Hamiltonian ([!]). 

The bulk Hamiltonian ^ possess a ground state which 
is four-fold degenerate. We label the different orthogonal 
ground states by a Greek letter, which can take the values 
^ = 0..3. This ground states are defined by the action of 
four operators T^, which act on the boundary of a state 
defined in the bulk in a similar way as ([6| . This state has 
spin 1 at each lattice site, but spin 1/2 at the boundary, 
then the action of the operators is to create spin 1 
also at the boundary. We have, for the Hamiltonian ([2| 



|GS^) = Tt(l,iV-l)|GS) 

N-l 

= Tt(l,7V-l)n(«Ifcli-al+ifcI)|0). (8) 

i=l 

The operators can be defined in term of the boson 
creation operators a and b as 



4(1, j) - aja] +blb], If (z,j) = albl+bla], 
rl(*,j)=*(al6]-6la]), T^it, j) = aU - bM. (9) 



then, we have i?Buik|GS^) = 0, (/i = 0,1,2,3). Linear 
combinations of this operators acting in the boundary of 
the chain can create separable states (take for example 
Tq + T3 ). Such states are also valid ground states of the 
chain. We will elaborate more on this point in section C 
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Negativity for the mixed system of 2 blocks 

We study the mixed system composed of two blocks A 
and B of length La and Lb, obtained by tracing away 
the lattice sites which do not belong to these blocks in 
the VBS ground state. This situation is described in Fig 

m 



B. Spin 1/2 at the boundary 



here we have introduced three types of tensors, the Kro- 
necker delta symbol in 4 dimensions S^, the diagonal 
tensor 5^1, — g^^" = diag(— 1, +1, +1, +1) and the Levi 
Civita tensor in four dimensions e^^pa, which is a totally 
antisymmetric tensor, with e^i,p^= sign of permutation 
(/i, ly, p, a) if (/i, V, p, a) is a permutation of (0, 1, 2, 3), and 
zero otherwise. 

Using this representation of the VBS state, it is easy 
to write down the density matrix ( 10 ) using the orthogo- 



nality of the bulk ground states, namely 



1 ... Lo-1 Lo ... m-1 m ... K K+1 ... M M + 1 ... N N + 1 

— O- ■ -O-^- ■ -G-^- ■ ■ ^>^- ■ ■ ■ ■ ■ ^3— 

'c (Lc^sites) A (U sites) D(L sites) B7usites)'^TTu'sites) 

FIG. 1: We made a partition of the VBS state in 5 sectors, 
labeled A,B,C,D and E as shown in the figure. To obtain 
the density matrix for the blocks A and B, we trace away the 
spin variables at the sites inside C, D and E. 



To define the blocks, we partition the TV + 2 sites of 
the chain into five different subsets, A,B,C,D and E, of 
different length. 

Given I,J,K, M, N five positive integers ordered as 
0<I<J<K<M<N+l,we define: 

• Block C = {sites i, < i < / —I}, with length 
Lc = I, 

• Block A = {sites i, I < i < J — 1}, with length 
La^J-I, 

• Block D = {sites i,J<i<K — 1}, with length 
L^K- J, 

• Block B = {sites i, K < i < M - 1} with length 
Lb = M - K and 

• Block E = {sites i, M < i < N + 1} with length 
Le = N + 2- M. 

We are interested in the density matrix for the mixed 
system of A and B blocks. We obtain this density matrix 
by tracing away the states on the C, D and E subspaces. 



{D^\D,) = Si^\^{L) with X^iL) = - + 



1 , z{L) 



(12) 



z{L)=l^--j ; 5^ = (-1,-1,3,-1), 

{C,E]T^{I~l,M)Tj{I~l,M)[C,E)=5ii. (13) 
we find that the density matrix pab is 



PAB — ^'^fii^pa 



(14) 



with the tensor M^^pcMauPa given explicitly by (sum- 
mation over dummy variables v and cr is assumed) 

+^z(L)5^"g"'5e^pA. [Spp - S^^) , (15) 

with S^a = (s^ + Sa)/2. This is the first explicit form 
of the reduced density matrix of two non complementary 
blocks for the AKLT model. This is important from the 
point of view of the fact that recently obtaining the re- 
duced density matrices of non-complementary blocks in 
a many-body system has been the focus of much interest 
[20] . We can identify two parts in ( 14 ) , the first term 



which does not depend on z and the rest which is linear 
in z. The first term is a projector on the ground states 
of the bulk of A and B, namely 



Po 



{A,B) = 5X\At.,Bp,){A^,B^ 



(16) 



while all the other terms, proportional to z{L), have van- 
ishing trace. If we call pi{A, B) to all the linear terms in 
z{L) on (14 1, we can write for brevity 



Pab = TrcDsip) 



(10) 



PAB = Po{A, B) + z{L)pi{A, B) 



(17) 



We can write the VBS state as a linear combination 
of products between the different four fold degenerate 
ground states of the bulk Hamiltonian ^ in the form: 
(implicit summation assumed) 



^ M^,p,Tl{I ~l,M)\C,A^,D,,Bp,E), with 



M^,p, = {-in6;gp,+d;g^ 



■ VOL 

ig e 



fiapa ) 



(11) 



From the expressions (14) and (15) we can obtain the 



partial transposed density matrix with respect to the A 
subsystem. 



6X + mK€~9'^g"^]s^.c. 



(18) 



+ iz{L)g^^' g"'^ eap\a{S pf) - Sf_,a) 



A^,Bp,){A^,Bp 
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If we perform a unitary transformation U on the basis 
vectors |Ap,i?p) defined by it's action on the basis as 
U\A^,Bp)=g^^-\A,,Bp), we find 



(19) 



iz{L)g^°'g''^ef^pXaiSpfi ~ S'^q) 



\A^,Bp,){A^,Bp\ 



from where, comparing with equations (14) and (15), we 
learn that 



Uplj,iz)U^ =PABi-z) 
With this resuh, we can state our main theorem: 



(20) 



Theorem 1. The negativity of the transposed density 
matrix p^g(z{L)) is strictly zero for two blocks separated 
by L>0. 

Proof. Consider the family of density matrices pab{z) — 
po{A, B) + z{L)pi[A, B), defined in eq. pT| Recalling 
that the space of density matrices is convex|T7] , meaning 
that for two density matrices pi, P2, the operator p — 
Xpi + (1 — X)p2 is also a density matrix for A € [0, 1], 
we proceed as follows. We take the first two members 
of the family pab{z), namely pab{zi) and pab{z2) for 
fixed zi = 1,Z2 = —1/3 [22] ■ By the convexity of the 
space of density matrices, p = \pab{zi) + (1 — X)pab{z2) 
is also a density matrix. Using ( |17[ ), we write explicitly 
p — pn{A, B) + (Azi + (1 — \)z2)pi{A, B). We can choose 
A = i(l - 3(-|)^) e [0, 1] for L>1. Using this A, we 
find 



P = Po{A,B)-z{L)pM.B). 



(21) 



Then p 



Pab{—z) is also a density matrix, for L >1. 
Now, by (20), p^si'^) ^1^° density matrix for z < 1 
(L > 0). 



"Een the negativity (sum of negative eigenval- 



ues) of p^g{z) vanish. 



□ 



For the case when the blocks are adjacent (L = 0), the 
negativity in the limit La, Lb ^ 1 is 



the ground state of a smaller chain. For this, first con- 
sider the 3 spin version of the above AKLT chain which 
has the ground state 



ialb\-a\bl){albl-albl)\0) 
n -1 1 + 2 + 



)|0)i 



U 



1, 



U 



l>i|- 2^2. 



(23) 



The 4x4 reduced density operator of the spins at sites 
and 2 (both being spin- 1/2 particles) is then 



1 



P02 



(I 



1 

2' 



1 



where 



+ 



2' 2'^ 2' 2' 
+ |^+)(V'+|),(24) 

)| -|- |)). It is easy 



v^d ' 2/1 2/ 

to check that po2 has a vanishing negativity. Now, the 
ground state of the 4 spin version of the AKLT chain with 
spin-l/2s at the ends may be obtained from the 3 spin 
version by bringing in a pair of spin- 1/2 particles 2' and 
3 in a singlet state \ip-) = + ^)\ ~ ^) - \ ~ ^)\ + ^)) 
and projecting 2 and 2' to their symmetric subspace to 
make a new spin-1 particle at site 2. After this operation, 
sites and 3 have spin-1/2 particles, while sites 1 and 2 
have spin-1 particles and the system as a whole is the 
ground state of the 4 spin AKLT model with spin-l/2s 
at the boundaries. The block composed of spins 2 and 3 
has thereby arisen from the spin 2 by a local action (a 
positive operator valued measurement) and a selection of 
a certain outcome (in this case a projection to a certain 
subspace). Since the negativity between and 2 in po2 
extracted from the 3 spin AKLT model is zero, then the 
negativity between and the block composed of spins 
2 and 3 in the 4 spin AKLT model should also vanish 
by virtue of the fact that negativity is an entanglement 
monotone (Eq.(17) of Ref. [21J). Continuing by induc- 
tion, adding spins on both sides of the chain to construct 
AKLT ground states with larger numbers of spins, it may 
thereby be shown that for blocks not touching each other 
the negativity vanishes. 



N{A,B)^^-^{z{LAf + z{LBr). (22) 

This means that as blocks get large, if they are touching, 
then their entanglement approaches that of a maximally 
entangled pair of q-bits. 

It is pertinent here to provide an alternative proof of 
the vanishing negativity of non-complementary blocks in 
the above AKLT ground state (i.e. with spin-l/2s at the 
ends) which may be obtained by exploiting the property 
that negativity is an entanglement monotone j21j , as well 
as the projective (filtering) mechanism through which the 
ground state of a larger AKLT chain can be grown from 



C. Spin 1 at the boundary 



The Hamiltonian in this case is just ([2| with no bound- 
ary terms. The ground state of a chain of spins 1, is four- 
fold degenerate as shown in 15 . We label these ground 
states by a greek index, which can take values from zero 
to three. We have 

N 

\VBSp) = Tl{l,N)Y[{a]bl, - al,bl)\Q). (25) 
1=1 



5 



1 2 ... m m+L-1 ... N-1 N 

* ^ ' 

L sites 

FIG. 2: Starting from one of the four ground states of the 
VBS chain, we trace away one block of L sites. This creates 
a mixed state of two blocks. 



The Tfj_ operators are defined as in ([9]). We define the 
states in this problem as A = 1 < sites < m — 1, B — 
m < sites < m + L — 1 and C = m + L < sites < N, 



\o)= n i<4bi+i-4+Am- 



(26) 



i,i+ieo 

The VBS state splits up in 



|VBS,) = Tlil,N)T^im-l,m)TUK,K + l)\A,B,C), 
= ^g^'-Tlil, N)Tl{m -1,K+ 1)| A, S, C.) (27) 

where if = m + L - 1, and \A,B,C) = \ A)\B)\C). 

The density matrix of the B subsystems can be com- 
puted directly from the expression (27), using the or- 



thogonality of the bulk states \Cu)- Defining |A, B)^, = 
TUm - 1,K + 1)\A,B), we have 



N = ^(3^9 - lOz(L) + 17z(L)2 + 5z(L) - 1), 
^ 1 + for L » 1. (31) 

In the thermodynamic limit La, Lb = oo (where the 
four ground states become indistinguishable [1 5j ) . the 
negativity quite surprisingly does not depend on the sep- 
aration of the blocks. The topology of the system de- 
termines the entanglement properties. Then when the 
blocks are adjacent (i.e L = 0), the negativity is maximal 
being A'^ = 3/2, while for any other separation N = 1/2. 

To obtain the previous result, we started from ( [27| , 
which is a entangled state between the first and the last 
spin. Taking linear combinations of |VBS) states, we can 
construct unentangled states between the first and the 
last spins. Relabelling ipj = Ui, ipf = bi, we can write a 
general unentangled state between the first and the last 
spin 1 /2 in the chain as 



r^VBS) ^ i;U'/Y[ialbl,- al4m. (32) 

i=l 

This states can be obtained as linear combination of 
the states |VBS^), for example, the state |^^VBS) — 
|VBSo) + IVBS3), and so on. The density matrix in this 
case, after tracing out block C is (c, d = 1, 2) 



PABit-i) = X^T^p{hN)\A,B)^{A,B\^Tpil,N). (28) 

We can write this density matrix in terms of the ground 
states of the A and B bulks. In terms of these we have 
(no implicit sum over /3) 



PabW) = Xa[(JxO-aCrp.Cri3][(Tx'CraCrf,'ap]\A^, Bx){Af,^ , By\, 

(Aq = Xa{L)), with [<J\aaCrf^ap] = Trace(crA0'aCr^cr/3) and 
5 being the complex conjugate of a. From this expression 
we can easily find the transposed density matrix respect 
to the A subsystem. We have 



PABic,d) = X^{L)4''^^/\A,B)^{A,Bl^i'i;N'' 



The partial transposed density matrix can be calculated 



A^(L)r'^Ai?)Mr'Ai?u- (33) 



directly from (34), the result is 



PAB{c,d)^u^{L)\-''A,B)^{'''A,Bl 



(34) 



withM^(L) = (A^(L)-i(-l)^(A2(£)-Ai(i)). Using this 
result, we compute the negativity of the system which 
turns out to vanish for L ^ 0, while for L = 0, we have 



PaB^P) = \a[(^\(yaCri_,ai3\[crx'<7a(T^i'(7l3\\A^_,f,Bx){A^,B 



(29) N{A,B-c,d) 



The negativity of the system is given by the sum of 
negative eigenvalues of ( 29 ) . Taking La, Lb and L the 



length of the blocks A, B and C respectively, we have in 
the limit of big blocks 



N ■ 



L=0 

1 

2 



3 

2 ' 



L>0 



for La,Li 



(30) 



liziLAf + ziLe)^) for La,Lb,L:^1. 



Another interesting limit is when one block is infinitely 
long compared with the other (La — 00, Lb = 1 for any 
separation L) 



y/l + z{LA)^ziLBr - z{LaY - z{LbY 
2~2^-dz{LA)z{LB) 

« i - \{z(La) + 0^"z(Lb))' {La,Lb » 1). 

(35) 

where we have introduced the symbol (jf^ which is equal 
to 1 if c = d and takes the value —1 if c 7^ d. 



D. Periodic boundary conditions 

Another case widely studied is the scenario with peri- 
odic boundary conditions. We can obtain this state from 
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the case studied in the previous section, antisymmetriz- 
ing the first and the last free spin 1/2 variables. This 
state is unique, as follows given that the coordination 
number for each spin is two II 65. 

In this state, we make a partition in four sectors, la- 
beled by their length as Li,La,L2,Lb, with Li + L2 + 
La + Lg = L the total length of the system. We trace 
away the states from the sectors that do not belong to 
AUB (See fig. ^. 




FIG. 3: We trace blocks 1 and 2, leaving a reduced density 
matrix in terms of the states of the A and B blocks. 



The reduced density matrix in this case is given by 

PAB^ M^fSa'(i'\A^,B^){A^,,Bfi,\. (36) 
The tensor Ma^a'p' is given explicitly by 

Mafia'p- = (5" 5^ Aa^(zi,Z2) ~~ ga^g" ^ aa' i~ Zi, - Z2) 
+ig^'^ g^" ei3\af,{Tal3a'0'{zi,Z2) - 5^ S'^,Taa'izi,Z2), 

(37) 

where zi = z{Li) and Z2 = 2(^2)- The tensors 
Aai3ix,y),Taa'{x,y) and Tq/jq'/j' (a;, y) are respectively 
given by 



1 + iSaSj3 + Sa+ Sp)xy 



1 + 3z(L) 

x + y 



and 



Tapa'f3'[x,y) = 4^i2z{L) ^^^y>- 

As with the case studied in section [Bj it is easy to 
see that if we write the partial density matrix in the 
form PAS (21, 2:2) = Pa + zipi + Z2P2 + 2122P3 the par- 
tial transposed density matrix is isomorphic (up to 
change of basis) to pab{—zi,—Z2)- We can use again 



the convexity of the space of density matrices to prove 
that the matrix pab{—zi,—Z2) is also a density matrix, 
for Li,L2 > 1. From this result we see that negativity 
vanish for Li,L2 > 1. 



Conclusions 

Understanding the entanglement structure of AKLT 
states is important because AKLT states have recently 
been recognized as useful resources for quantum informa- 
tion processing [181 US] ■ We have computed the entangle- 
ment between non-complementary blocks of spins in such 
a system. We have found that the AKLT chain allows 
an analytic computation of the negativity between non- 
complementary blocks in a short range (realistic) spin 
model, which has not been accomplished in the other 
models studied to date. 

The entanglement properties of a system of two blocks 
made up from the ground state of an AKLT system de- 
pends dramatically on whether the blocks touch each 
other. We have proved that the negativity of a system 
of two blocks separated by L > 1 sites obtained from a 
unique ground state ([6]), is strictly zero. 

In the case when all the sites contain a spin 1 and 
open boundary conditions, the entanglement properties, 
and particularly the negativity, depends on the config- 
uration between the sites at the boundary. Each of the 
spin 1 variables can be constructed from the symmetriza- 
tion of two spin 1/2 states. In particular, the sites at 
the boundary can be constructed through this process. 
In this boundary sites, however, one of those spin 1/2 
can be in any state, i.e. | f) or | |). We can construct 
Bell states between the first and the last spins 1/2 in 
the chain. This states have non vanishing negativity re- 
gardless the separation between the subsystems A and 
B inside the chain. If, on the other hand, the spins 1/2 
at the end of the chain are in one of the four separa- 
ble states (i.e. calling the first and the last site 1 and 
N respectively, | titiv), I tiiiv),! Iitw),! IiIn)) then 
the negativity vanish when the subsystems A and B are 
separated. Thus, in practice, as there is no easy way to 
isolate one of those AKLT ground states which have a 
Bell state between its end spins, we would expect vanish- 
ing negativity even for the case when all sites contain a 
spin 1. 

We have also shown that the vanishing negativity be- 
tween separated blocks of spins also holds for the periodic 
AKLT which has a unique ground state. Other impor- 
tant results of our work include a general expression for 
the reduced density matrix of an arbitrary pair of blocks 
of spins in a AKLT chain that includes the case of sep- 
arated blocks, and the negativities between those non- 
complementary blocks that touch each other at a point. 
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